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ABSTRACT 



The commonality of collisionally replenished debris around main sequence 
stars suggests that minor bodies are frequent around Sun-like stars. Whether or 
not debris disks in general are accompanied by planets is yet unknown, but debris 
disks with large inner cavities - perhaps dynamically cleared - are considered 
to be prime candidates for hosting large-separation massive giant planets. We 
present here a high-contrast VLT/NACO angular differential imaging survey 
for eight such cold debris disks. We investigated the presence of massive giant 
planets in the range of orbital radii where the inner edge of the dust debris 
is expected. Our observations are sensitive to planets and brown dwarfs with 
masses >3 to 7 Jupiter mass, depending on the age and distance of the target 
star. Our observations did not identify any planet candidates. We compare 
the derived planet mass upper limits to the minimum planet mass required to 
dynamically clear the inner disks. While we cannot exclude that single giant 
planets are responsible for clearing out the inner debris disks, our observations 
constrain the parameter space available for such planets. The non-detection of 
massive planets in these evacuated debris disks further reinforces the notion that 
the giant planet population is confined to the inner disk (<15 AU). 



Subject headings: circumstellar matter - planetary systems - stars: individual (HD 
105, HD 377, HD 107146, HD 202917, HD 209253, HD 35850, HD 70573, HD 25457) 
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Introduction 



CoUisionally re plenished de 
Sun-like stars (e.g. iMeyer et al 



)ris d ust surrounds about 10-20% of the main sequence 



20071 ). Such widespread evidence for minor body collisions 



demonstrates that planetesimals orbit most stars. It is natural to ask whether or not rocky 
and giant planets are also present in these systems. No convincing co rrelation could yet 



be fo und bet ween close-in exoplane ts and the presence of debris (e.g. 



Moro-Martfn et al. 



2007. but see 



Beichman et al. 



20051 ). However, the presence of massive giant planets has 



been often invoked to account for the observed azimuthal or radial 
radii in many debris disks (e.g. 



Greaves et al. 



offers several alternative mechanisms (e.g. 



2005 



Wilner et al 



asym metries at large 



Takeuchi Sz Artymowicz 



2002h. While theory 



2001 



Wvatt 



200,^1 ) 



dynamical clearing of dust parent bodies by giant planets remains a feasible and excitin g 



theoretical possibility (e.g. 



Moro-Martfn fc Malhotra 



2007 



2005 



Morbidelli et al. 



20071 ). 



Quillen 



2006 



Levison et al. 



Examples for such possibly dynamically-clea red disks inc . 



disks around the you ng Sun-like stars HD 105 flMeyer et al 



ude two recently identified 



( iWilliams et al. 



20041 ) and HD 107146 



20041 ). Both disks were found to exhibit strong excess emission at 
wavelengths longer than 30 /im, while displaying no measurable excesses shortward of 
20 /im. The detailed analysis of the spectral energy distribution of HD 105 suggests that it 
is consistent with a narrow du st ring (<4 AU) w ith an inner radius of ^ 42 AU, if the du st 



grains emit like black bodies (IMeyer et al. 



20041 ). Using a similar model 



Williams et al. 



(120041 ) showed that the excess emission from HD 107146 is consistent with arising from 
cold dust (T=51 K) emitting as a single-temperature black body. The lack of measurable 
infrared excess shortward of 25 fim illustrates that the inner disk regions are well cleared 



of dust: for HD 107146 t 
(T=51 K 



lereis at most 140 x less warm dust (T=100 K) than cold dust 



Williams et al. 



20041 ). The findings of the spectral energy distribution model for 
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HD 107146 have been confirmed by direct imaging with the Hubble Space Telescope, that 
strengthen the cas e for a large featureless dust ring outside of an evacuated inner cavity 



( lArdila et al. 



2004h . 



Recent high-contrast imaging su rveys have hinted on the general scarcity of giant 



planets at such large separat ions (e.g. 



2007 



Lafreniere et al 



Masciadri et al. 



2005 



Kasper et al. 



2007 



Biller et al. 



20071 ). Quantitative statistical analysis of the non-detections 



demonstrates that - at a 90% confidence level - the giant planet population cannot extend 
beyond 30 AU if it follows a r°-^ radial distribution, consistent with the radial velocity 
surveys. Th e statistical analysi s suggests an outer cut-off for the giant planet population 



at <15 AU (IKasper et al. 



20071 ). If so, dynamically cleared cold debris disks may be the 



ssignposts for rare large-separation giant planets, ideally suited for direct imaging studies. 

In this paper we report on a VLT/NACO high-contrast imaging survey for large- 
separation giant planets around HD 105, HD 107146, and six other similar disks. In 
the following we will review the target stars and disks, the observations, followed by a 
comparison of our non-detections to lower planet mass limits set by dynamical clearing 
simulations. 



1.1. Targets 



Our targets were selected from the sample of 328 Sun-like stars (0.7-2.2 Mq) targeted 
in the Formation and Evolution of Planetary Systems Spitzer Space Telescope Legacy 



program (FEPS, 



Meyer et al. 



20061 ) . From this sample we identified 8 southern stars, 
which: a) display strong infrared excess emission at long wavelengths (A > 20/im); b) no 
measurable excess emission at shorter wavelengths; and, C) are young and close enough to 



permit the detection of planetary-mass objects within the inner radius of the cold debris. 
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Table [T] gives an overview of the key parameters of the target stars. The typical lower mass 



limit for the debris in the systerns is 



of our Kuiper-belt (IMeyer et al. 



to 10 Mm, making these disks massive analogs 



2004: iKim et al. 



2005 



Hillenbrand et al., in prep). The 



disks of HD 105 and HD 107146 — included in our sample — have inner evacuated region s 



2004 



Williams et al. 



2004J ). 



with an estimated radii of ~ 40 AU and ~ 31 AU (IMeyer et al.l 
The other six disks exhibit spectral energy distributions similar to HD 105 and HD 107146. 
Based on the similarity of the excess emissions and the almost identical spectral types all 
eight disks are expected to have cleared-out inner disks of similar size. The only possible 
exception in this sample is HD 202917, for which the re-calibration of the IRAC fluxes after 
our VLT observations revealed a faint, but likely real infrared excess even at wavelengths 
shortward of lO/im, suggesting that this inner disk may harbor small, but non-negligible 
amounts of warm dust. 

In the following we discuss briefly the results of the age determination for these sources 
as this has direct impact on the sensitivity of our observations to giant planets. A more 
detailed discussion of the ages of the whole FEPS sample will be presented in Hillenbrand 
et al. (in prep). We briefly summarize the upper and lower age estimates (tmm and tmax) 
for each star along with the most likely age tprob, where available. HD 105 has already 
reached the main sequence {tmin='^7 Myr) and its chromospheric activity suggests a tmax of 
225 Myr (Hillenbrand et al., in prep.). Ve ry likely a member of t he Tuc-Hor moving group 



(IMamajek et al 



20041 ) its tprob is 30 Myr flHoUenbach et al 



20051 ). HD 377 is also a main 



sequence star (tmm > 25 Myr) and the chromospheric activity suggests that tmax =220 Myr. 
The median of four other age indicators sets tprof)=90 Myr. For HD 107146 we adopt the 
age range of 80-200 Myr. HD 202917 is a likely member of the Tuc-Hor moving group 
{tmin = tprob = 30 Myr) and its upper age limit is set by its Li-abundance, higher than 
that of the Pleiades (tm n,^ < 100 Myr). H D 35850 is suggested to be a /5 Pic Moving Group 



member (tmm=12 Myr, 



Song et al 



20031 ) and its observed rotation rate sets a reliable 
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upper age limit of tmax=^00 Myr (Hillenbrand et al. in prep.; cf. 



Barnes 



20071 ). HD 70573 



is among the few star s that are known to harbor both a debris disk and a giant planet. 



Setiawan et al 



(120071 ) found an m2sini = 6.1 Mj possible planet on a 1.76-AU orbit. 



A combination of different age indica tors suggest a tmm = 30 Myr for HD 70573 and a 



tprob = 60 Myr; 



Setiawan et al 



(120071 ) quotes tmaa; = 125 Myr. Based on Li-abundance and 
chromospheric activity, position on the color-magnitude diagram and the analysis of its 
space motions Mamajek et al. (in prep.) estimates that HD 209253 has t^m = 200 Myr and 



tmax < 1-6 Gyr with t 



prob 



member of the A 


3 Dor moving group 


Zuckerman et al. 


2004 


)• 


Luhman et al. 



giving a very strong lower age limit (tmin = 50 Myr 
( 120051 ) derives an age of 75-125 Myr (we adopt tprob = 75 Myr), while the upper age limit is 
set by the chromospheric activity (tmax = 170 Myr). 



2. Observations and Data Reduction 



Our 8 targets were observed with ESO's Very Large 



adaptive optics system (ILenzen et al. 



2003 



Rousset et al. 



Teles cope using the NACO 



20031 ). The observations were 



carried out in service mode in late 2006 and early 2007. The weather conditions were 
excellent with typical visual seeing of O'.'S and clear skies. 

We used the spectral differential imaging mode (SDI) of NAC O in order to enhan ce the 



Lenzen et al 



200J). The 



contrast for any methane-rich cold (T< 1200 K) companion (e.g. 
SDI mode uses two Wollaston prisms to split the incoming light rays into four beams of 
nearly identical light path. These rays pass through four narrow-band filters, two of which 
are identical. The three different filters (/i, /2 and /s corresponding to 1.575, 1.600, and 
1.625 /im) probe the 1.62 /im methane feature and the adjacent continuum. Because the 
SDI mode uses the 1024 x 1024-pixel S13 camera of NACO with a 13" x 13" field of view, 
the simultaneous acquisition of 4 images in this field reduces the effective field of view to 
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about 3" X 3" 



The achieved contrast, however, is not as good as predicted, probably due to the 
combined effect of read-out noise and a slightly lower Strehl ratio. With the reduced 
contrast our observations were sensitive only to planets beyond the 1-3 Mj planet mass 
range. These planets - at the young ages of out targets - are too hot to display the 1.62 



/xm-methane feature (IBurrows et al 



20031 ). The lack of the methane feature results in 



almost identical planet fluxes in the /i, /2 and /s filters, rendering the SDI technique 
inefficient. 

Instead, we opted to reduce the data taken in the /i filter (Ac = 1.575/im, 6\ = 0.025/im) 
in the angular diff erential imaging mode, i.e. without applyin g the spectral differ ential 



imaging step (e.g.. 



Mueller fc Weieelt 



1987 



Marois et al. 



2006 



Kasper et al. 



data r educ tion was performed with a dedicated pipeline, as described in detail in 



fl2005h and 



Janson et al 



2007 ). The 



Kellner 



(120071 ). The frames taken at a given rotator angle were averaged 
and the collapsed frame corresponding to one angle was subtracted from the other. This 
procedure cancels out residual static or quasi-static features from the instrument, whereas 
any companion will remain as a combination of a positive and negative point source. The 
intensity of the residuals (at a certain separation from the primary) is characterized by 
taking the standard deviation in a 9 x 9-pixel square (O'.'llxO'.'ll ~ (2.7A/D)^ ) centered 
on that separation, at 180 evenly sampled angles, and taking the median of the results. 
This is repeated for all separations to create a radial profile of the error distribution in 
the image. When combined with the brightness of the prim ary, this yields the achieved 



Janson et al. 



2007 for more 



contrast as a function of separation in the final image (see, 
details). Previous artificial planet tests on identical data sets processed with the reduction 
pipeline used here sh owed that 3-sigma- bright sources would have been reliably identified 



as candidate planets (jJanson et al. 



20071 ). Lacking any such detection, we used 3-sigma 
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fluxes as upper limits on the brightness of any companions to the target stars. In order 
to convert the achieved contrast in fi to the more commonly used H-band, we derived a 
conversion factor by comparing the flux densities in the two filters in a si mulated spectru m 



of a giant planet in the age and mass range probed by our observations (IBurrows et al. 
20081 ) ■ 



3. Results 

The NACO/ADI observations acquired high-contrast, high-resolution images of the 
8 target stars and their immediate environment (typically between 10 and 70 AU). In spite 
of the sensitive observations we could not identify giant planet candidates or any other 
point sources in the images. 

Our data analysis allows us to set firm upper limits to the brightness of the sources 
that would have been identif ied as a candidate. We use the age estimates in Table [T] and 



planetary evolution models (iBaraffe et al. 



20031 ) to convert the achieved sensitivities to 



planet masses. These limits are shown in Fig. [T] as a function of separation from the target 
stars for the lower and upper age limits of our stars. 

Using the high-contrast observations we can study the probable range of radii in these 
disks at which the dust debris and the parent body planetesimals reside. Our measurements 
exclude the presence of any brown dwarf companions for virtually all our targets at these 
radii. For HD 209253, the oldest star in our sample, the images exclude companions down 
to the brown dwarf/giant planet mass boundary (13 Mj ) at radii 20 AU or greater. For 
the the youngest source HD 35850 (12-100 Myr) our observations exclude any giant planet 
companions down to 3 -4 Mj betw e en 25 to 45 AU, if the star belongs to the (3 Pic moving 



group as suggested by 



Song et al. 



(120031 ). For the other six sources our observations are 
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Tablc 1. Target parameters. 



Target 


R.A. (J2000) 


Dec. (J2000) 


V-mag." 


Dist. [pc]"^ 


Sp. Type 




HD 105 


00 05 52.6 


-41 45 11 


7.51 


40 


GOV 


27 Myr / 30 Myr / 225 Myr 


HD 377 


00 08 25.7 


+06 37 01 


7.59 


40 


G2V 


25 Myr / 90 Myr / 220 Myr 


HD 25457 


04 02 36.8 


-00 16 08 


5.38 


19 


F5V 


50 Myr / 75 Myr / 170 Myr 


HD 35850 


05 27 04.8 


-11 54 03 


6.30 


27 


F7/8V 


12 Myr / 12 Myr / 100 Myr 


HD 70573 


08 22 50.0 


+01 51 34 


8.69 


70 


G1/2V 


30 Myr / 60 Myr / 125 Myr 


HD 107146 


12 19 06.5 


+16 32 54 


7.04 


29 


G2V 


80 Myr / - / 200 Myr 


HD 202917 


21 20 50.0 


-53 02 03 


8.65 


46 


G5V 


30 Myr / 30 Myr / 100 Myr 


HD 209253 


22 02 33.0 


-32 08 02 


6.63 


30 


F6/7V 


200 Myr / 500 Myr / 1.6 Gyr 



^AU magnitudes and distances from the Hipparcos catalog, except for the distance of HD 70573, which is a main 
sequence-distance. 



'The age estimates are discussed in the text. 



Table 2. Log of the observations. 



Target 


UT Dates 


NDIT X DIT'^ 


Frames per Angle 


Angles 


On-Sour ce 


Strehl 


HD 105 


07/20/06+08/18/06+08/21/06 


1 X 4s 


144 


0°, 33° 


19 min 


43% 


HD 377 


08/27/06 


1 X 4s 


144 


0°, 33° 


19 min 


37% 


HD 25457 


08/12/06+08/13/06 


1 X 4s 


144 


0°, 33° 


19 min 


46% 


HD 35850 


08/13/06 


1 X 4s 


144 


0°, 33° 


19 min 


40% 


HD 70573 


03/02/07 


21 X 5s 


16 


0°, 33° 


56 min 


58% 


HD 107146 


04/260/06+0 05/26/06 


1 X 4s 


95 


0°, 33° 


13 min 


51% 


HD 202917 


06/23/06 


1 X 12s 


32 


0°, 33° 


13 min 


45% 


HD 209253 


07/09/06+07/16/06+07/23/06 


1 X 4s 


144 


0°, 33° 


19 min 


49% 



^DIT — Detector integration time; NDIT — number of integration averaged on— chip. 
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typically sensitive to ~6 Mj at orbital radii >30-40 AU. 

Note, that the dominant uncertainty of the upper limits stems from the difficulty of 
stellar age determination and from the poorly constrained initial conditions for giant planet 



accre tion phase, giant planets may start with mu ch lower luminosities (e.g. 



20071 ) than assumed by the hot-start models (e.g. 



Baraffe et al 



2003 



i (e.g. 


Mar lev et al. 


Burrows et al. 


2003) 



4. Discussion: Inside— Evacuated Debris Disks Without Massive Giant 

Planets? 

If single planets are responsible for clearing out the inner disks in the observed systems, 
they will be located very close to the inner edge of the debris. Although the available data 
does not allow the direct measurement of the inner disk radii, simple black body fits to 
the spectral energy distributions provide reliable lower limits. We adopt these limits from 
Hillenbrand et al. (in prep.) and note that they range from 6.2 to 28 AU. Fig. [1] shows the 
limits for the individual sources (dashed lines). We stress again that these are lower limits 
— the real inner disk radii are probably somewhat larger. 

Given the range of detectable planet masses in Fig. [H we assess whether dynamical 
clearing by a less massive and therefore undetectable planet could still be a feasible 
mechanism to explain the lack of measurable quantities of warm dust, or on the contrary, 
if th e dynamical clearing scenario may be rejected. Becaus e our disks are devoid of 



gas ( 



'ascucci et al. 



Moro-Martm fc Malhotra 



20061 ) we use the dynamical models by 
( l2005l ) to investigate the effect of a giant planet on the dust population. These models 
investigate the efficiency of dust particle ejection by gravitational scattering as a function 
of planet mass and planet location. In these models the dust particles are released from 
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Discovery Space for the Survey Targets 




20 40 20 40 60 

Distance from the Star [AU] 



Fig. 1. — The discovery space of the NACO observations (shaded), the minimum predicted 
planet masses required for a single planet to scatter out > 90% of the planetesimals at 
a given radius (dotted lines) and the lower limits for the disk inner radius (dashed lines, 
Hillenbrand et al., in prep.). The presence of planets within the shaded parameter range 
is excluded by our observations. The upper and lower shaded sensitivity curves mark the 
limits for a possible younger and older stellar age. For each figure the outer radius of the 
field of view is given in astronomical units in the lower right-hand corner. Note, that the 
reversing sensitivity curve for HD 70573 and HD 209253 is due t o the reversal of th e planet 



mass-luminosity curve in the corresponding evolutionary phase (iBaraffe et al 



20031). 
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an outer belt of planetesimals and drift inward toward the central star under the effect of 
Poynting-Robertson drag, scattering as they cross the orbit of the planet and naturally 
creating a dust-depleted region inside its orbit. However, as shown in the upcoming study 
by Hillenbrand e t al. (in prep .) all our disks are collision-dominated, as are bright debris 



disks in general (IWyatt 



20051 ). This means that the dust particles may not have time to 



drift too far from the parent bodies before getting eroded by collisions down to the blow-out 
size. Thus, in these collision-dominated disks the dust generally traces the location of 
the planetesimals. Therefore, we need to evaluate the effect of a giant planet on the 



planetesimal population rather than on the dust 



particles. Because gravitationa. 



Moro-Martm fc Malhotral (120051 ) are also 



scattering 



is a process independent of mass, the models of 
applicable to planetesimals as long as these can be considered to be "test particles" (i.e. 
their masses are negligible with respect to that of the planet). 

Using the above models we evaluate what is the mass of the least massive planet 
that can open a gap in the planetesimal distribution. In order to provide a good model 
for the evacuated inner gaps we require that the planet scatters out at least 90% of the 
planetesimals. A dotted line in Fig. [1] shows these lower masses as a function of orbital 
radius. We find that for planets in the 5-30 AU range, a planet mass of at least 2-5 
Mj masses is required. Thus, it is conceivable that giant planets in the mass range 2 to 
5 Mj clear the gaps and still remain undetected by our survey. However, we point out that 
in the cases of most of our targets, and in particular for HD 35850, the parameter space 
that such a planet can occupy is limited. 

Given that our high-contrast imaging survey did not find single, large-separation giant 
planets that m ay be respon si ble fo r clearing the inner disks, we briefly explore alternative 



mechanisms. iBesla fc Wd ( 120071 ) provides a useful summary of the proposed models 



and we will only highlight here a particularly interesting proposal proposed recently by 
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Moro-Martin et al.l (120071 ) for the system HD 38529, where two close-in planets trigger 
secular resonances that affect the planetesimal population in the outer disk. In this model 
the eccentricity of the planetesimals at the location of the secular resonances is excited, 
thus enhancing the rate of collisions and truncating the planetesimal disk. 

Given the range of possible mechanisms that may lead to the formation of dust rings it 
is probable that in-depth studies of the individual systems will be required for judging the 
feasibility of the proposed models on a case-by-case basis. However, our non-detections 
show that the presence of evacuated large inner holes in cold disks could be related to 
localized dust production or concentration of dust grains by dust-gas interactions rather 
than dynamical clearing by single massive giant planets. 



5. Conclusions 

We present results from a high-contrast angular differential imaging survey of 8 cold 
debris disks, selected to have significantly or totally evacuated inner disks. Our observations 
searched for massive giant planets that may be responsible for carving out the inner holes 
in the observed cold debris disks. For most of our targets we reach typical sensitivities of 3 
to 7 Mj between 20 to 50 AU separations, but did not identify any likely planet or brown 
dwarf candidates. 

By comparing the derived planet mass upper limits to lower limits derived from 
dynamical scattering models (typically 2-5 Mj between 10 and 30 AU), we limit the 
parameter space available for any single planet capable of efficiently clearing out the inner 
planetesimal disks. 

Our survey complements recent direct imaging surveys of nearby young stars indicating 
that massive giant planets at large separations are very rare. Cool debris disks with large 
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inner evacuated cavities remained promising possible exceptions to this rule until now. 
However, the combination of our observational upper limits and theoretical lower limits 
strongly suggest that massive giant planets at large separations are not present in most of 



these systems, reinforcing the finding that the outer cut-of f for the giant 



is probably at 15 AU or at even smaller semi-major axes (iKasper et al. 



plan et distribution 



2003) 
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the Office of Space Science to the Life and Planets Astrobiology Center (LAPLACE). We 
would like to thank members of the FEPS team for their help in characterizing the target 
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